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Abstract 

We construct a matrix model that reproduces the topological string partition 
function on arbitrary toric Calabi-Yau 3-folds. This demonstrates, in accord with 
the BKMP "remodeling the B-model" conjecture, that Gromov-Witten invariants of 
any toric Calabi-Yau 3- fold can be computed in terms of the spectral invariants of a 
spectral curve. Moreover, it proves that the generating function of Gromov-Witten 
invariants is a tau function for an integrable hierarchy. In a follow-up paper, we will 
explicitly construct the spectral curve of our matrix model and argue that it equals 
the mirror curve of the toric Calabi-Yau manifold. 
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1 Introduction 

In the topological string A-model, the object of study is the moduli space of maps from 
a Riemann surface £ 9 of genus g to a given Calabi-Yau target space X. Its partition 
function is the generating function of Gromov-Witten invariants of X, which roughly 
speaking count these maps. 

In recent years, deep connections have been unrooted between the topological string on 
various geometries and random matrix models. A classic result in the field is that inter- 
section numbers, which are related to the Gromov-Witten theory of a point, are computed 
by the Kontsevich matrix integral pQ, see also [2]. In the Dijkgraaf-Vafa conjecture [3] 
such a connection is obtained between the topological B-model on certain non-compact 
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Calabi-Yau manifolds and a 1-matrix model. A novel type of matrix model [I] inspired 
by Chern-Simons theory is associated to the topological string in [5], yielding matrix 
model descriptions of target spaces obtained from the cotangent space of lens spaces via 
geometric transition. This work is extended to chains of lens spaces and their duals in 

In the 20 years that have passed since topological string theory was formulated [TJ [8] , var- 
ious techniques have been developed for computing the corresponding partition function. 
The topological vertex method [H] solves this problem completely for toric Calabi-Yau 
3-folds at large radius, furnishing the answer as a combinatorial sum over partitions. 
On geometries with unit first Betti number (the conifold and 0{— 2) — > CP 1 x C), this 
formalism yields the partition function as a sum over a single partition with Plancherel 
measure. In [ID] , such a sum was rewritten as a 1-matrix integral. More complicated 
examples, such as the topological string on geometries underlying Seiberg-Witten SU (n) 
theory, can be written as sums over multiple partitions (TTJ [12} [13] . 1-matrix integrals that 
reproduce the corresponding partition functions were formulated in [H] . Multi-matrix in- 
tegrals have arisen in rewriting the framed vertex as a chain of matrices integral [15] . Its 
Hurwitz-numbers limit (infinite framing of the framed vertex geometry) was shown to be 
reproduced by a 1-matrix model with an external field in [TBI E] • 

Here, generalizing the method of [10] , we are able to formulate a matrix model which 
reproduces the topological string partition function on a certain fiducial geometry, which 
we introduce in the next section. Flop transitions and limits in the Kahler cone relate 
the fiducial geometry to an arbitrary toric Calabi-Yau manifold. As we can follow the 
effect of both of these operations on the topological string partition function, our matrix 
model provides a description for the topological string on an arbitrary toric Calabi-Yau 
manifold. 

By providing a matrix model realization, we are able to transcribe deep structural insights 
into matrix models to the topological string setting. E.g., our matrix model involves a 
chain of matrices, and chain of matrices integrals are always tau functions for an inte- 
grable system. Our matrix model realization hence proves integrability of the generating 
function of Gromov-Witten invariants. Moreover, matrix models satisfy loop equations, 
which are known to be equivalent to W-algebra constraints. A general formal solution 
to these equations was found in p2], centered around the introduction of an auxiliary 
Riemann surface, referred to as the spectral curve of the system. The partition and corre- 
lation functions of the matrix model are identified with so-called symplectic invariants of 
this curve [TH]. The BKMP conjecture [20], building on work of (21], identifies the spec- 
tral invariants of the mirror curve to a toric Calabi-Yau manifold with the topological 
string partition function with the Calabi-Yau manifold as target space. In a forthcoming 
publication [22], we will compute the spectral curve of our matrix model explicitly, thus 
establishing the validity of this conjecture. 

Finally, we would like to emphasize that many different matrix models can yield the same 
partition function (justifying the choice of indefinite article in the title of this paper). 
An interesting open problem consists in identifying invariants of such equivalent matrix 
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models. A promising candidate for such an invariant is the symplectic class of the matrix 
model spectral curve. 

The outline of this paper is as follows. In section [2j after a very brief review of toric 
geometry basics, we introduce the fiducial geometry and the notation that we will use in 
discussing it throughout the paper. We also review the transformation properties of the 
topological string partition function under flop transitions, which will relate the fiducial 
to an arbitrary toric geometry, in this section. We recall the topological vertex formalism 
and its application to geometries on a strip [23] in section [3] Section [4] contains our main 
result: we introduce a chain of matrices matrix model and demonstrate that it reproduces 
the topological string partition function on the fiducial geometry. By the argument above, 
we thus obtain a matrix model description for the topological string on an arbitrary toric 
Calabi-Yau manifold, in the large radius limit. We discuss implications of this result in 
section |5j and point towards avenues for future work in section [6} 



2 The fiducial geometry and flop transitions 

Toric geometries present a rich class of very computable examples for many questions in 
algebraic geometry. The topological vertex formalism provides an algorithm for computing 
the generating function for Gromov-Witten invariants on toric 3 dimensional Calabi-Yau 
manifolds. These are necessarily non-compact and have rigid complex structure. 

The geometry of toric manifolds of complex dimension d can be encoded in terms of a d 
dimensional fan S, consisting of cones of dimensions to d. We denote the set of all n 
dimensional cones as E(n). Each such ra-cone represents the closure of a (C*) d ~ n orbit. 
In particular, 1-cones correspond to hyper surf aces, and for d = 3, our case of interest, 
2-cones correspond to curves. 

The fan for the class of geometries we are interested in is constructed by triangulating 
a finite connected region of the Z 2 lattice containing the origin, embedding this lattice 
in Z 3 within the (x,y) plane at z = 1, and defining the cones of the fan via half-lines 
emanating at the origin and passing through the vertices of this triangulationfj] 

We can associate a dual diagram to such toric fans, a so-called web diagram, spanned by 
lines orthogonal to the projection of 2-cones onto the Z 2 lattice. In the web diagram, the 
relation between the dimension of the components of the diagram and the submanifold of 
the toric geometry they represent coincide: 3-cones (points) correspond to vertices, and 
2-cones (curves) to lines, see figure HI 

1 The canonical class of a toric manifold is given by the sum over all torically invariant divisors. The 
construction sketched above guarantees that this sum is principal, hence the canonical class trivial: the 
monomial associated to the 1-cone (0,0, 1) generates the class in question. See e.g. [23] . 



3 



Figure 1: Example of a box triangulation, corresponding to a 3 dimensional toric fan. The diagram 
in red is the dual web diagram. Vertices of the triangulations (faces of the web diagram) correspond to 
1-cones, edges correspond to 2-cones, and faces (vertices of the dual) correspond to 3-cones. 



2.1 The fiducial geometry 

The geometry X we will take as the starting point of our considerations is depicted in 
figure [2] 
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Figure 2: Fiducial geometry Xo with boxes numbered. 

Since the torically invariant curves play a central role in our considerations, we introduce 
a labeling scheme for these in figure [3] (i, j) enumerates the boxes as in figure § and we 
will explain the a-parameters further below. 




Figure 3: Labeling curve classes, and introducing a-parameters. 

In the following, we will, when convenient, use the same notation for a torically invariant 
curve E, its homology class [E] e H2(Xq, Z), and its volume or associated Kahler parame- 
ter J s J, given a Kahler form J on X . The classes of the curves r^j, s^-, introduced in 
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figure |3| are not independent. To determine the relations among these, we follow [25, page 
39, 40]. Consider the integer lattice A spanned by formal generators e p , with p e S(l) 
1-cones of the toric fan, 

a = {£ 

pes(i) 



A p e p |A p e Z} 



Each torically invariant curve, corresponding to a 2-cone of the fan, maps to a relation 
between 1-cones, and thus to an element of the lattice A, as follows: a 2-cone a is spanned 
by two integral generators v\ and t>2, and it is contained in precisely two 3-cones, which 
are each spanned by Vi, v% and one additional generator t> 3 , t> 4 respectively. These vectors 
satisfy the relation Ylt=i ^i v i = 0' where the Aj can be chosen as relatively prime integers, 



and as v 3 and v 4 lie on opposite sides of a, we can assume that A 3 , A 4 > 0. [22] shows 
that on a smooth variety, the sublattice A^ generated by the elements J2t=i ^« e « °f A is 
isomorphic to H2(Xq,7j). We call this isomorphism A, 



X:H 2 (X ,Z)^A h . 



Figure [4] exemplifies this map. 




vi : (0,1,1) 
v 2 : (1,0,1) 
^3 : (0,0,1) 
«4 : (1,1,1) 



A = (-1,-1,1,1) 



Figure 4: The 2-cone a corresponds to the relation A among 1-cones. 

It allows us to easily work out the relation between the various curve classes. Consider 
figure [5] 



«6 



<'3 



V 7 







\ v 4 


\ti+l,j-l 

ri,j-i\ 



Vl 



V 2 



Figure 5: Determining the relation between curve classes. 
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The images of the curve classes depicted there under A are, 



A(r M ) 


= e 5 + e 6 


- e 4 


- e 7 


Hrij-i) 


= e 2 + e 3 


- ei 


- e 4 




= e 3 + e 7 


- e 4 


- e 6 




= ei + e 5 


- e 2 


- e 4 



We read off the relation 



By symmetry, we also have 



(2-1) 



+ s i+l,j-l ■ 



A moment's thought convinces us that this constitutes a complete basis for the space of 
relations. We can solve these in terms of the classes of the curves r i; Sj, tij, i,j = 0,l,... 
depicted in figure [6j which hence generate H 2 (Xq, Z). The explicit relations are 




mi ■ • ■ 1 L n,m 



^0,1)^1,1) ■ ■ i tn,\ 
^0,0)^1,0) • ■ • j ^n,0 

fo n r 2 r 3 r n 

Figure 6: Fiducial geometry with choice of basis of H 2 (X , Z). 



fc=l 



fc=l 



Our computation for the partition function on 3£o will proceed by first considering the hor- 
izontal strips in the toric fan describing the geometry, as depicted in figure [2j individually, 
and then applying a gluing algorithm to obtain the final result. 

For each strip, we find it convenient to write the curve class wjj e H 2 (3io, Z) of the curve 
extending between two 3-cones which we label by I and J (recall that 3-cones correspond 
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to vertices in the dual web diagram), with J to the right of J, as the difference between 
two parameters a/ and aj associated to each 3-cone, 



wu = ai-aj. (2.2) 

We call these parameters, somewhat prosaically, a-parameters. It is possible to label the 

curve classes in this way due to their additivity along a strip. In terms of the notation 
introduced in figure [3j we obtain 



By invoking the relation (2.1 ), we easily verify that upon gluing two strips, the curve class 
of a curve extending between two 3-cones / and J on the lower strip is equal to the class 
of the curve between the 3-cones I' and J' on the upper strip, where the cones / and I' 
are glued together, as are the cones J and J', 

Wu = wpj> . (2.3) 

This allows us to identify the parameters aj = ap and aj = aj> associated to 3-cones 
glued together across strips. 

Note that the basic curve classes s$ are not captured by the parameters a it j. 



2.2 Flop invariance of toric Gromov-Witten invariants 

Under the proper identification of curve classes, Gromov-Witten invariants (at least on 
toric manifolds) are invariant under flops. Assume X and X + are related via a flop 
transition, : X — > X + . In a neighborhood of the flopped (—1,-1) curve, the respective 
toric diagrams are depicted in figure [7j 



T2 r 2 



















T"3 




n t 3 





Tl 



T 4 T 4 



Figure 7: X and X + in the vicinity of the (-1,-1) curve. 

The 1-cones of S^, corresponding to the toric invariant divisors of X, are not affected by 
the flop, hence can be canonically identified with those of X + . The 2-cones Tj in these 
diagrams correspond to toric invariant 2-cycles C i: Cf in the geometry. The curve classes 
of X push forward to classes in X + via 

^([C ]) = -[C + ], ^([a]) = [C+] + [C + ]. (2.4) 
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All other curve classes of X are mapped to their canonical counterparts in X + . Under 
appropriate analytic continuation and up to a phase factor (hence the oc in the following 
formula), the following identity then holds [261 E31 127] . 

Zgw\%i Qoj Qii ■ ■ ■ i Qii Q) Zgw(% > 1/Qoj QaQii ■ ■ ■ i QoQ&i Q) , (2-5) 

i.e. 

GW g (X, Q ,Qi,..., Q 4 , Q) = GW g (X + , l/Q , Q Q U Q Q Q±, Q) . 

Any toric Calabi-Yau manifold X with Kahler moduli Q can be obtained from a sufficiently 
large fiducial geometry (X , Q ) upon performing a series of flop transitions and taking 
unwanted Kahler moduli of Xq to oo. Once we obtain a matrix model reproducing the 
topological string partition function on the fiducial geometry, extending the result to 
arbitrary toric Calabi-Yau 3-folds will therefore be immediate. 

As an example, we show how to obtain the P 2 geometry from the fiducial geometry with 
2x2 boxes in figure [8} 




Figure 8: We obtain local P 2 from the fiducial geometry with 2x2 boxes by performing five flops and 
then sending the Kahler parameters of the unwanted edges to oo. 



3 The partition function via the topological vertex 
3.1 Gromov-Witten invariants 

Gromov-Witten invariants A/" 9i d(3C) roughly speaking count the number of maps from a 
Riemann surface of genus g into the target space X, with image in a given homology class 
D = (D 1: . . . , D k ) G H 2 (X, Z). They can be assembled into a generating series 

GW g (X,Q) = J2KAX)Q D - 

D 
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Each GW g (X,Q) is a formal series in powers Q D = YliQi ' °f the parameters Q = 
{Qii Q2, ■ ■ ■ , Qk)i the exponentials of the Kahler parameters. 

We can introduce a generating function for Gromov-Witten invariants of all genera by 
introducing a formal parameter g s (the string coupling constant) and writing 

00 

GW(X,Q,g s ) = 9 2 s 9 ~ 2 GW g (X,Q). 

g=0 

It is in fact more convenient to introduce disconnected Gromov-Witten invariants J\f* D (X), 
for possibly disconnected surfaces, of total Euler characteristics x, and to define 

z 6 w(x,q,9s) = e cw ^ = j2q d K, D m- 

D X 

For toric Calabi-Yau manifolds, an explicit algorithm was presented in [§] for computing 
Z GW via the so-called topological vertex formalism, proved in [2S1 EH]- 

3.2 The topological vertex 

In the topological vertex formalism, each vertex of the web diagram contributes a factor 
C q (a,fl, 7) to the generating function of GW-invariants, where the a, /3, 7 are Young 
tableaux associated to each leg of the vertex, and C 9 (a,/3,7) is a formal power series in 
the variable q, where 

q = e~ 9s . 

Topological vertices are glued along edges (with possible framing factors, see [9j) carrying 
the same Young tableaux a by performing a sum over a, weighted by , with Q encoding 
the curve class of this connecting line, 

Young tableaux« e edgese verticesu=(ei,e2,e3) 

Note that in practical computations, the sum over representations can ordinarily not be 
performed analytically. A cutoff on the sum corresponds to a cutoff on the degree of the 
maps being counted. 

The equality 

Zgw(%, Q,9s) •^vertex(-^j Qi Q) 

holds at the level of formal power series in the Q's, referred to as the large radius expansion. 
It was proved in |29J that the log of the right hand side indeed has a power series expansion 
in powers of g s . 
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3.3 Notations for partitions and q-numbers 



Before going further in the description of the topological vertex formula, we pause to fix 
some notations and introduce special functions that we will need in the following. 



3.3.1 Representations and partitions 



Representations of the symmetric group are labelled by Young tableaux, or Ferrer dia- 
grams. For a representation 7, we introduce the following notation: 

• 7»: number of boxes in the i-th row of the Young tableau associated to the repre- 
sentation 7, 71 > 72 > • • • > 7^ > 0. 

• The weight I7I = Ylili'- the total number of boxes in the corresponding Young 
tableau. 

• The length ^(7): the number of non- vanishing rows in the Young tableau, i.e. 7i = 
iffi>Z(7). 

• The Casimir ^(7) = J7 liili — 2z + 1). 

• 7 T denotes the conjugate representation, which is obtained by exchanging the rows 
and columns of the associated Young tableau. We have |7 T | = I7I, l( r f T ) = 71, and 
K (7 T ) = -«(t)- 

An integer d > will denote a cut-off on the length of representations summed over, 

1(7) < d- 

Most expressions we are going to write will in fact be independent of d, and we shall argue 
in [22], following the same logic as in [10J based on the arctic circle property [30], that 
our results depend on d only non-perturbatively. 



To each representation 7, we shall associate a parameter a as introduced in (2.2). 

Instead of dealing with a partition 7, characterized by the condition 71 > 72 > 
0, it will prove convenient to define the quantities 

Kil) =Ji-i + d + a, 

which satisfy instead 

hi > h 2 > h 3 > ■ ■ ■ > h d > a . 
The relation between 7 and ^1(7), for the off-set a = 0, is depicted in figure [9| 
We finally introduce the functions 

ar.( 7 ) = g M7) . 

In terms of the ^(7), we have 



>id> 

(3.1) 



where C d>a = \{d - l)(2d - 1) + a(a + 2d - 1). 
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h h h h h 

10 9 8 7 6 



3 2 1 



Figure 9: Relation between a partition 7 and h^). 

3.3.2 q-numbers 

We choose a string coupling constant g s such that the quantum parameter q = e~ 
satisfies \q\ < 1. A g-number [x] is defined as 



[x] — q 2 — q 2 = 2 sinh 



xg s 



(3.2) 



g-numbers are a natural deformation away from the integers; in the limit q — >■ 1, — [x] — > 
x. 

We also define the g-product 



00 1 



n=l 



- -q 

X 



The function g(x) is related to the quantum Pochhammer symbol, g(x) = [q/x; q]^, and 
to the g-deformed gamma function via T q (x) = (1 — qY~ x g{l) / gi^q 1 '*) . g(x) satisfies the 
functional relation 



g(qx) 



x 



For T q , this implies T q (x + 1) = T q (x), the quantum deformation of the functional 
equation T(x + 1) = xT(x) of the gamma function, which is recovered in the classical limit 
q — > 1. The central property of g(x) for our purposes is that it vanishes on integer powers 
of q, 

g(q n ) = if n G N*. 
Moreover, it has the following small In q behavior, 



H n=0 



■D n B r , 



111 



(\nq) n Li 2 _ n (l/x) 
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where Li n (x) = X^kLi * s ^ ne polylogarithm, and B n are the Bernouilli numbers 

Bq = 1 , Bi = —- , B 2 



1 

6 



B2k+i — if k > 1 (see the appendix). 

We shall also need the following function f(x) 



1 



/(*) 



g(x)g(q/x) -fain. 

p Jlnq p In q 

s(l) 2 

— In g 



t7T H-K *> 

In (jr ' In g ^ 



lnx 1 



27T — 2«7T 



lng 2 lng' lng 



where # is the Riemann theta-function for the torus of modulus —2in/ lng. This relation- 
ship is the quantum deformation of the classical gamma function identity 

e~™ x /r(l - x)T(x) = sin (irx)/iT . 



3.4 The partition function via the vertex 

We begin by considering a single horizontal strip of the fiducial geometry, as depicted in 



figure 10 



Pi Pi 



T 



a. 




Pi 



a 



T 

n 



0o 



01 



Figure 10: A horizontal strip of the fiducial geometry and its corresponding web diagram. 

Of the three legs of the vertex, two point in the direction of the strip and connect the 
vertex to its neighbors. One leg points out of the strip, either above or below. This leg 



carries a free representation, or /?• in the notation of figure 10 The partition function 



will hence depend on representations, one per vertex (i.e. face of the triangulation). 
A note on notation: since each 3-cone carries a representation (which up to the final 
paragraph of this subsection is held fixed) and an a-parameter (see figure [3]), we will 
identify the a-parameters by the corresponding representations when convenient. 
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Using the topological vertex, it was shown in [23] that the A-model topological string 
partition function of the strip is given by a product of terms, with the individual factors 
depending on the external representations and all possible pairings of these. Applied to 
the fiducial strip, the results there specialize to 

[A > «i ] Q ft lli<y [«< > Z 3 / ] Qctf,/3j \Pi > «j ] Q ft , Qj 

We explain each factor in turn. 

• Each vertex 7 = a, or 7 = /3 4 T contributes a representation dependent factor to the 
partition function, which we have denoted by [7]. It is the n — > 00 limit of the Schur 
polynomial evaluated for x% = q?~\ i = 1, . . . ,n, given explicitly by 



n (^-^) n 

l<j<j<d <=1 



1/(1) 

= A(X( 7 )) e-^f^^W-^) e ~wi(x(« i 

We recall that ^(7) = 7 i — z + d + a 7 , and we have defined = q hi and the diagonal 
matrix X(~f) = diag(q 1 ,q 2 ,...,q d ). Furthermore, A(X) denotes the Vandermonde 
determinant of the matrix X, 

A(X)= J] (xj-Xi), 

l<i<j<d 



and we have written 

U(X,a) = -g s ln 



9i£ 
90) 



(\nX) 2 

Ui(X,a) = K ' - (a + d- l)lnX \nq + C(a,d), 

aja+d-T) , (<l-l)(2d-l) 



where C(a, d) , r> 

We have 



[7] = ] b T ] , i^{i T ) = > 

[ 7 T ] = A(X( 7 )) e"^ trC7(x(7) ' a ^ e-57 tr ^ (x(7) ' a ^, 



and thus 



where 



C/"i(X,o) = ^lnX lnq + C(a,d). 
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C a ^ is another constant which depends only on a and d and which will play no role for 
our purposes. 

• In addition, each pair of representations contributes a factor, reflecting the contribution 
of the curve extended between the respective vertices. In the nomenclature of |23j, the 
representations cti are all of same type, and of opposite type relative to the If we take 
i < j, representations of same type (corresponding to (-2,0) curves) contribute a factor 
of 

[ai, aj] or /?,•] , 

whereas representations of different type (corresponding to (-1,-1) curves) contribute a 
factor of 

1 1 

or 



The pairing is given by [3H H21 1321 123] 



d 7i j d bi oo 

n n [a a5 + ^ _ < + 5j n n [a as _ J + % _ rf] ru^v 

i=l j=l L 7 J i=l j=l 1 T J fc=0 



= '-^n* ri(-^ 

i=i l ' «=i 

11 y -j--'- g(q a i~ as ) g(q as ~ a t) 

OC A(X( 7 ),X(5)) e -^ Tr ^(7)^) e -^TrC/(X(5),a 7 ) g -i (Tr U 2 (X( 1 ),a s )+Tr C/ 2 (V(5),a 7 )) ^ 

(3.4) 



where the square brackets on the RHS denote g-numbers as defined in (3.2), the symbol 
A(X(j),X(8)) signifies 

A(X( 7 ),X(5)) = IJ(^-%)) = II(^ (<5) -^ (7) )> (3-5) 

and 

U 2 (X,a)=Q, 

(\ n X) 2 1 

f/ 2 (X,a) = ^ '- - (a + d- -) lnX lng + ivrlnX 

2 2 
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The parameter Q lt § reflects, given a choice of Kahler class J of the metric on 3£ , the 
curve class of the curve C extended between the vertices labeled by 7 and 5 via 

w 1>5 = / J , Q ltS = q Wl ' s ■ 



Jc 

By the definition of the a-parameters, 



Substituting these expressions into (3.3), we obtain 

Zstrip(ttcb ■ ■■ i a n', 0q i ■ ■ ■ , ) = 

lI,, / A(.V(n,)..V( 
x TTe~^ tr(ys(x(ai)) ~^ (Jf(ai))) TTe"^ trV ' l(x(Q!)) 



Y[ e ~^ tT{Vsixiai)) ~ v t {xiai))) JJ. 

i i 



where we have denoted by a = (a , a 1; . . . , a n ) (resp. b = (b , b\, ... , &„)) the a-parameters 
of representations on the upper side (resp. lower side) of the strip, and defined 

n 

V s (X) = -g s J2^(9(q aj /X)), (3.7) 

j=0 

and 



V t (X)=lnXlnq ( 1 - £ 

( £ ~ Z^' ~ a i) ) 



K(-X") = lnX lng 



3.5 Gluing strips 

To obtain the partition function for the full multistrip fiducial geometry Xq, we must glue 
these strips along the curves labelled Sjj in figure [3j 
Denoting the representations atj^ on line % collectively by 

this yields 

m+l n m 

^vcrtcx(Xo) = Z (n , m) (a m+ i, al)= II Z strip(ai, c^) JJ JJ fel . 

ctj i, j=0,...,n; i=l,...,m «=1 j'=0 i=l 

(3. 

Our goal now is to find a matrix integral which evaluates to this sum. 
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4 The matrix model 



4.1 Definition 

Consider the fiducial geometry j£ of size (n + 1) x (m + 1), with Kahler parameters 
ti,j = a it j — a ij+ i, rij = dij+i — (k+ij, and Sjj, as depicted in figures [3] and [6j We 
write 

0>i — (Oo,i) • • • , 

Assume that the external representations are fixed to a m+ i = (ao, m +i, a i,m+i, ■ ■ ■ , «n,m+i) 
on the upper line, and c? — ( a o,o> a i,o, ■ ■ ■ > a n,o) on the lower line (for most applications, 
one prefers to choose these to be trivial). 

We now define the following matrix integral 2mm (mm for Matrix Model), 

m+l „ m+l 



^MM(Q,^,« m +i,«o) = A(X(a m+1 )) A(X(a )) TT / dM* TT / di^ 

i=o J^MiVi) i=1 JH n (R+) 

m m 



i=l i=l 
m+l m 



i=l i=l 



)tr In Mj 



e trln/ (Mo) gtr ln/ m+ i(M m+1 ) f[p trh /i( M i) 



■TTgtrln/iCMi) _ 

i=l 

All matrices are taken of size 

N=(n+l)d, 

where d is the cut-off discussed in section 13.3.11 We have introduced the notation 

X(a m+1 ) = diag(X(a m+1 ) i ) i= i,... )iV , X(a m+1 ) jd+k = gM«i,m+i) , 

X{a ) = dmg(X(a )i)i=i,...,N , X(a ) jd+k = q h ^°\ 
for k — 1, . . . , d, j — 0, . . . , n. A(X) = Yli<j(Xi — Xj) is the Vandermonde determinant. 



Vg^x) was introduced in (3.7). For % — 1, . . . , m, we have defined 



/<(*) = n 



._ ( . ^(ig 1 a ^)g(q a ^/x) 

The denominator of these functions induces simple poles at x = q a ^ i+l for j = 0, . . . , n and 
/ G Z. The numerator is chosen such that they satisfy the relation fi(qx) = fi(x). This 



16 



enforces a simple / dependence of the residues taken at x — q a i' i+ \ given by a prefactor q l 
- a fact which will be important in the following. These residues are in fact given by 



0(1)2 e (|+&)(l+«J.i-«*,i)lng — 
lco^j,^ y « 11 g(^qaj,i-ak^ (I — q a k,i- a 3,i"jg(^q a k,i-aj,i"j 



,-,.-) ( - -I-,'--. ! [ I-!-,; .-„, . ! hi,; — - 

Res /,(./•) 7" I -'I"--' TT 

(4.2) 

where is independent of the integer /. 
The parameters Si are defined by 

<Si = So,j-l + to,i-l — s j,i-l ~ + • (4-3) 

fc<i k<j 



The final equality holds for arbitrary j, and can be verified upon invoking (2.3) repeat- 
edly. 

For i = and i — m + 1, we define 

/o(z) 



n;=onii(*-^ o) ) 
1 

/m+l(^) 



n;= nli(z-^- +i) ) 

Notice that if the representations «o or ct m+ i are trivial, i.e. hi(atj t o) = d — i + cij t0 or 
^i(«i,m+i) = d — i + aj )Tn+ i, we have 

f M _TT g( g g 1 ~ aj, °" d ) f - TT g(gg 1 " aj ' ,m+1 ~ d ) 

respectively. The functions /o and f m +i have simple poles at x = q h i( a i,o) (resp. x = 
g^iCaj.m+i)) for Z = 1, . . . , d, with residue 

/j, ;i = Res / (x 



^ko) nti(? ft!(aj ' o) - g hi( ^'' o) ) n^(^ i(a - o) - ^ ( °- o) ) 

i i 



fj,m+l;l — R es fm+l( x 



The / dependence here is more intricate than above, but this will not play any role since 
the partitions a^o and <x, im +i are kept fixed, not summed upon. 

The integration domains for the matrices Ri are Hn(M.+ ), i.e. the set of hermitian ma- 
trices having only positive eigenvalues. For the matrices Mj, % = 1, . . . , m, the integration 
domains are H^iTi), where 

n 

i n • 

3=0 
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7j,i is defined as a contour which encloses all points of the form g a J> i+N ; and does not 
intersect any contours jf.,u Cm) ^ (k,l). For this to be possible, we must require that 
the differences — ayy be non- integer. The normalized logarithms of two such contours 



are depicted in figure 11 




Figure 11: Two contours surrounding points a + N and b + N, such that a — b^TL. 
We have defined 

H N {T i ) = {M = UAU\ UeU(N), A = diag(Ai, . . . , X N ) GT,}, 

i.e. H^(Ti) is the set of normal matrices with eigenvalues on Tj. By definition, the 
measure on H^^Ti) is (see [33J) 



dM = — A(A) 2 dUdA 



(4.4) 



where dll is the Haar measure on U(N), (normalized not to 1, but to a value depending 



only on N, such that the Itzykson-Zuber integral evaluates as given in (4.6) with pre- 



factor 1), and dA is the product of the measures for each eigenvalue along its integration 
path. 

The integration domains for the matrices M , M m+ i are Hn(Tq), HN(T m+ i) respectively, 
where 

n n 
3=0 j=0 



(4.5) 



The goal of the rest of this section is to prove that the matrix integral (4.1) reproduces 



the topological string partition function for target space the fiducial geometry j£ . 



4.2 Diagonalization 

Let us first diagonalize all matrices. We write 

Mi = UiXiU}, 



R i = U i Y i U, l 



t 

i J 



where Ui and Ui are unitary matrices. 
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By the definition (4.4), the measures dMi and dRi are given by 



dMi = ^ A(Xi) 2 dUi dX t , dR i = ^ A{Y i fdtj i dY i . 
The matrix integral thus becomes 

zMQ,g s ,a m+1 ,^) = A(x( 7';, ) L A 3 (x( " o)) n / ^A^n / 

l ■) j =0 J^i i=1 Jr^ 

m+l m+1 

Hdu, 

i=0 i=l 

m m 

J] e ^ tr JJ e ^ tr [v 3i _ 1 (x i _ 1 )-y s .(x,_ l) ] 

8=1 1=1 

n+ 1 m 
JJ ^trXiUlUiYiUjUi e ^\*X l - 1 vt_ 1 U i Y i & i V i - 1 JT ( 



i=l i=l 
m+l 

'-trXiUlUiYiUlUi p =^tr X^ul^U^ulU,. , TT e (5l+ ^ )tr A 

i=l i=l 

m 

e tr ln/o(X ) gtr ln/ m+ i(X m+ i) TT g tr ln/i(X 4 ) _ 



n 

i=l 

Next, we introduce the matrices C7j, for i = 1, . . . , m + 1, via 

We can express {To, ... , £7 m +i, and Ui, . . . , U m +i, in terms of these matrices and U m+ i, 

U% = C/ m +i £/ m +i C/ m +i C/ m L/ m . . . Z7j+i t/j+i , 

C^i = f^m+l f^m+1 f^m+1 • • • "i+l ■ 

With this change of variables, we arrive at 

iWQ,&,<Wi,<£) = A(X( ^n ) Lts (X( " 0)) II / dX 4 A(X 4 ) a n / ^A^) 5 

/m+l m+l 
rff/m+i n ^ n ^ 

i=l i=l 
m m 

JJ e ^ tr [^ (Xl) - ys -i (Xl) ] TT e ^ tr K_i(^-i)-^(^-i)] 



i=l i=l 

m+l m 



i=l i=l 

m 

e trln/ (X ) g tr ln/ m+ i(X m+ i) TT g tr ln/i(Xj) 



n 

i=i 

Notice that the integral over C/ m+ i decouples, and J dU m+ i = Vol(U(N)). 
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4.3 Itzykson-Zuber integral and Cauchy determinants 

The JJi and t/j appear in the form of Itzykson-Zuber integrals 



'(•">/ jt,£tro ^S' (4 - 6) 

where x p and y q are the eigenvalues of X and Y. We thus have 

zmi(Q,g,,cL+i,c%) oc A(x( 7 i ;;jL+3 (x((?o)) n / ^wn'u,^ 

m m 

t=l i=l 
m+1 _. m 

n a-^,^) /( — Xi-uYi) n> ( 



3 (5, + ^)tr lnX t 



,=i 9s 9s 



e trln/ (X ) g tr ln/ m+1 (X m+1 ) TT g tr In /, (X, ) 



i=l 

m+1 » m+1 



OC 



n 

A(-Y(a„, +1 )) A(X(3 )) W , -p, , 

i^ 3 Mir,"- M V 

m 

A(X )A(X m+1 ) gei^W-^W] 

i=l 

m m 

i=l i=l 
m+1 

FT det(e^ (x ^ (y ^) det(e^ (x '-^ (y) <) 



e (S*+^)tr InX 



, P,9 P,9 
i=l 

m 

e trln/ (X ) e trln/ m+1 (X m+1 ) TT g tr In fi{Xi) 



i=l 



where we have dropped an overall sign, powers of g s , and the group volume Vol(U(N)) 
which are constant prefactors of no interest to us. 

Next, we perform the integrals over Yj along 



dYdet(e^^ h(Yh ) det(e 

N 



ly poo 
a a p=l J ° 

£E(-'H-') 8 n (I . lh y ' (x 

a 5 p=l ^*-1^(p) 



i)a(p) 
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Nig? det 



1 



p,g \(Xi-i) p - (Xi) q 

Note that the integral is only convergent for (X i ) -( p ) — (Xj_i)g.( p ) < 0. For Xj that violate 
this inequality, we will define the integral via its analytic continuation given in the third 
line. 

An application of the Cauchy determinant formula, 



det 

yields 



1 \ X\l<i<i<n( X 3 ~ X i)(Vj ~ Vi) 



x i + Vj J l<i<j< n Yli,j=l( X i + Vj) 



A(Xj_x, Xi 



where the notation A(Xj_i,Xj) was introduced in (3.5). Evaluating the Yi integrals thus, 



and continuing to drop overall signs and powers of g s , our matrix integral becomes 

7 in - A(X(dWi)) A(X(q )) V? f 2 
Zm&{Q,g 8 ,oi m+ i,a ) oc (W) m + 3 H / ' " 

i=0 



L 

L A(X,_ 1 ,X, 



i=l i=l 

m+1 m 

TT - TT (S l + ^)trlnX I 

11 A(Y. , 11 



i=l 

m 



, tr In/uiAul , t, ln./„. . , ( ..V . , ,) g tr In /» (X, ) 



= 1 



4.4 Recovering the sum over partitions 

Following the steps introduced in [Hj in reverse, we next decompose the diagonal matrix 
Xi into blocks, 

Xi = diag (Xo,i, Xij, . . . , X nj j) , 

where each matrix Xj i is a d x ci diagonal matrix whose eigenvalues are integrated on the 
contours 7^ surrounding points of the form q a ^ i+N . We arrive at 



~ (n . -»T\ A(X(q m+1 )) A(X(a )) TtA / .„ 
Zmm{Q, g s ,a m+ i,a ) oc (N\) m + 3 1111/ ^ 



8=0 j=0 " (Tj.*)" 

m+1 



AfXWX 1 n A(X,_ 1 )A(X i 
A(X )A(X m+1 ) 11 A(Xj_i,Xj) 
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Tfl 1(1 

1=1 

m m 



i=l i=l 

m m 

e trln/o(X ) e trln/ ro+ i(X m+1 ) ]~|" g tr In "pj gOSi+fj) tr In 



t=l i=l 

with 

A(X i _ 1 )A(X i ) _ MX,, :jlI f Ai.V,,) n^Afe-O^ j I ; , ; A( A ,,. .V,, ) 
A(X i _ 1 ,X i ) n^A(X^_ l5 X M ) 

Our next step is to evaluate the dXjj integrals via Cauchy's residue theorem. The poles 
of the integrands lie at the poles of fa, and the zeros of A(Xj_i,Xj). However, we have 
been careful to define our contours 7^ in a way that only the poles of fa contribute. 
These lie at the points g a ^ +N . Hence, the integrals evaluate to a sum of residues over the 
points 

{X hi ) t = q"^^)> , 
where each (hj^)i is a positive integer. 

Since the integrand contains a Vandermonde of the eigenvalues of Xj^, the residues vanish 
whenever two eigenvalues are at the same pole of fa, i.e. if two (hj } i)i coincide. Moreover, 
since the integrand is symmetric in the eigenvalues, upon multiplication by N\, we can 
assume that the (hj^i are ordered, 

(Mi > (hih > (hj,ih > ■ > (hj,i) d > 0. 
The (hj t i)i hence encode a partition via (hj^i = (aj^i—i+d, and we have reduced our 



integrals to a sum over partitions. In terms of the function hi(a) introduced in (3.1 ), 

= g h '(«j,0 , hi(a jti ) = {hj ti )i + a jti , 
hi(atj,i) > h 2 (a jti ) > ■ ■ ■ > h d (a jti ) > a^. 

Notice that unlike fa, % = 1, . . . , m, fo and / TO +i only have a finite number of N = (n+ l)d 
poles. Since the (hj t0 )i, (hj im +i)i respectively can be chosen pairwise distinct and ordered, 
/o and fm+i act as delta functions in the integrals over the N x N matrices X and X m+ x, 
and fix these to the prescribed values X(c? ) and X(a m+ i) respectively. 

Performing the integrals hence yields 

Z MM (Q,g s ,d m+ i,dl) oc A(X(a m+1 )) 2 A(X(a )) 2 

v ^ A(X(^_ 1 ))A(X(^)) 

{ aj ,i\j=0,...,n;i=l,...,m+l} i=l V V i; ' V 1/1 
m m 

JJ tr K(X(a l ))-V Si _ 1 TT g ^ tr [Vfc^ (X(3i_i))-Va, (X(<Si-i))] 

8=1 8=1 
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Notice that 



Furthermore, 



m m+1 n d 

ii nun i 

i=i 



1717 Res / 

1111 h,(a, ) 



\ n h l( a j,i) 

i=0 j=0 2=1 V 

1 



, ^Ko) A(X(a )) 2 

1 



n n ^ 



2 <? 



m+1 



A(X(a m+1 )) 2 



Res ./; = q h ^ /,, , 



where /^j computed in (4.2) is independent of hi(atj 7 i). We thus have 



n d 



g&trtaXfSO TJJJ Res J . =e (S i +l)trlnX(5 i ) TT 
J= 2=1 V"^' / j=0 



Upon substituting the expression (4.3) for Si, we finally arrive at 



oc 



(4.7) 



-2^MM(<3, 9s, «m+l, <5q ) 
m n 

mi e 

i=l jf=0 {a!j j i|j=0,...,n;i=l,...,m+l} 



m+1 

n 

i=l 
m 



n, ApTfe-O) n, A(X(a il4 )) n,-<i A(X(a i)4 _i), X(a M -i)) A(X(a,,), X(a M )) 



i_l)-g)tr lnJf(a M ) 



i=l i=l fc=0 

m m n 

JJ e ^ tr K-i^^-i))-^^^-!))] JJ JJ e (5-Ei<*(tti,i-aj,*4-i))te i^K,.) 

i=l i=l fc=0 

m n 



nrr 

i=l i=0 



Comparing to (3.6) and (3.8), we conclude 



m+1 



^mm(Q ; 5^; <3 m +i, <3q ) oc 2^ W Z 

ctj i,j=0,...,n;i=l,...,m i=l 



i=o i=l 



i.e. 



Z MM (Q,g s ,a m+1 ,a^) oc Z vcrtcx (X ) = e E 9 ^- 2 Gty 9 (* ) 
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Up to a trivial proportionality constant, we have thus succeeded in rewriting the topolog- 
ical string partition function on the fiducial geometry X as a chain of matrices matrix 
integral. 



By our reasoning in section 2.2 



this result extends immediately to arbitrary 
toric Calabi-Yau 3-folds as follows. We have argued that any such 3-fold can be obtained 
from a sufficiently large choice of fiducial geometry via flops and limits. The respective 
partition functions are related via (2.5). Upon the appropriate variable identification, we 
hence arrive at a matrix model representation of the topological string on an arbitrary 
toric Calabi-Yau 3- fold. 



5 Implications of our result 

We have rewritten the topological string partition function as a matrix integral. This 
allows us to bring the rich theory underlying the structure of matrix models to bear on 
the study of topological string. 

The type of matrix integral we have found to underlie the topological string on toric 
Calabi-Yau 3-folds is a so-called chain of matrices. This class of models has been studied 
extensively [351 E2L an d many structural results pertaining to it are known. 

5.1 Loop equations and Virasoro constraints 

The loop equations of matrix models provide a set of relations among correlation functions. 
They are Schwinger-Dyson equations; they follow from the invariance of the matrix inte- 
gral under a change of integration variables, or by an integration by parts argument. 

Loop equations for a general chain of matrices have been much studied in the litera- 
ture, in particular in [361 EH EHJ [39]. They can be viewed as W-algebra constraints (a 
generalization of Virasoro constraints) |40j . Having expressed the topological string parti- 
tion function as a matrix integral, we can hence conclude that Gromov-Witten invariants 
satisfy W-algebra constraints. 

Moreover, a general formal solution of the loop equations for a chain of matrices matrix 
model was found in [39] , and expressed in terms of so-called symplectic invariants F g of 
a spectral curve. The spectral curve for a matrix integral is related to the expectation 
value of the resolvent of the first matrix in the chain, 

/ 1 \ (0) 
W(x) = (tr — ) . 

V ; \ x-M / 

The superscript indicates that the expectation value is evaluated to planar order in a 
Feynman graph expansion. The symplectic invariants F g (C) of an arbitrary spectral curve 
C were defined in [19]. [39] proved that for any chain of matrices integral Z, one has 

g 
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with C the spectral curve associated to the matrix integral. 



Calculating the spectral curve of a chain of matrices matrix model with complicated 
potentials poses some technical challenges. We will present the spectral curve for our 
matrix model (4.1) in a forthcoming publication [22J. 



5.2 Mirror symmetry and the BKMP conjecture 

The mirror I of a toric Calabi-Yau 3-fold X is a conic bundle over C* x C*. The fiber 
is singular over a curve, which we will refer to as the mirror curve S% of X. It is a plane 
curve described by an equation 

<% : H(e x , e y ) = 0, 

where if is a polynomial whose coefficients follow from the toric data of X and the Kahler 
parameters of the geometry. 

Mirror symmetry is the statement that the topological A-model partition function with 
target space X is equal to the topological B-model partition function with target space 
X. 

Extending work of Marino |21j proposing a relation between the formalism of [19J and 
open and closed topological string amplitudes, Bouchard, Klemm, Marino and Pasquetti 
(BKMP) conjecture in [20] that 

GW g (X) = F g {S k ) . 

Here, the Fg's are the symplectic invariants introduced in [19]. The main interest of this 
conjecture is that it provides a systematic method for computing the topological string 
partition function, genus by genus, away from the large radius limit, and without having 
to solve differential equations. 

This conjecture was motivated by the fact that symplectic invariants have many intriguing 
properties reminiscent of the topological string free energies. They are invariant under 
transformations 5—7-5 which conserve the symplectic form dx A dy = dx A dy, whence 
their name [TH]. They satisfy holomorphic anomaly equations [41 J, they have an integrable 
structure similar to Givental's formulae [421 H3J HU HS1 HE], they satisfy some special 
geometry relations, WDVV relations J4T] . and they give the Witten-Kontsevich theory as 
a special case [19j HH] . 

BKMP succesfully checked their claim for various examples to low genus. 
The conjecture was proved for arbitrary genus in [10] for X a Hirzebruch rank 2 bundle 
over P 1 (this includes the conifold). Marshakov and Nekrasov [13] proved Fq = GWq for 
the family of SU(n) Seiberg-Witten models. Klemm and Sulkowski [14] . generalizing [TU] 
to Nekrasov's sums over partitions for SU(n) Seiberg-Witten gauge theories, proved the 
relation for Fq, building on work in [49]. In fact, it appears straightforward to extend 
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their computation to arbitrary genus F g . In [SO], Sulkowski provided a matrix model 
realization of SU(n) gauge theory with a massive adjoint hypermultiplet, again using a 
generalization of [TO] for more general sums over partitions. Bouchard and Marino [51] 
noticed that an infinite framing limit of the BKMP conjecture for the framed vertex 
X = C 3 implies another conjecture for the computation of Hurwitz numbers, namely that 
the Hurwitz numbers of genus g are the symplectic invariants of genus g for the Lambert 
spectral curve e x = ye~ y . That conjecture was proved recently by another generalization 
of [TO] using a matrix model for summing over partitions [IE] , and also by a direct cut and 
join combinatorial method [52J. The BKMP conjecture was also proved for the framed 
vertex X = C 3 in [53J [51], using the ELSV formula and a cut and join combinatorial 
approach. 

Since we have demonstrated that the topological string partition function is reproduced 
by a matrix model, we can conclude that the Gromov-Witten invariants coincide with the 
symplectic invariants 

9 G 

with C the spectral curve of our matrix model. We will compute C explicitly in a forthcom- 
ing work |22j, and demonstrate that it indeed coincides, up to symplectic transformations, 
with the mirror curve S%, thus proving the BKMP conjecture for arbitrary toric Calabi- 
Yau 3-folds, in the large radius limit. 

5.3 Simplifying the matrix model 

The matrix models associated to the conifold or to geometries underlying Seiberg-Witten 
theory have a remarkable property: the spectral curve is the same (perturbatively and up 
to symplectic transformations) as the one of a simpler matrix model with all ^-functions 
replaced by only the leading term in their small In q expansion. We will demonstrate in a 



forthcoming work [22] that this property also holds for our matrix integral (4.1). We can 
hence simplify the potentials of our matrix model, arriving at 



m+1 „ m+1 



Z Bimp (Q,g s ,a m+ i,6%) = A(X(a m+1 )) A(X(a )) TT / <IM, ] [ / dlU 

i= o JBHiFi) i=1 JH n (n+) 



W e9s 



f T,]= O (Li 2 (q^-yM l )-Li 2 (q a 0,i-^/M i )) 



i=l 
m—1 



n e3s 



tr £" =0 ("2(g a ^ /M 1 )-Li 2 (g a ^ l+1 /Mi)) 



i=0 

m+1 m 
i=l i=l 



where the matrix Mj is of size rii = ^ ■ rijj . 
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Classical limit 



In the classical limit, the dilogarithm Li 2 becomes the function x\nx, and we have 

m+1 „ m+1 „ 

Z eS . cl (Q,g s ,a m+1} c%) = A(X(a m+1 )) A(X(a )) TT / r/.U, j | / dRi 

i=0 JH^CCi) i=1 JHn(R+) 



n 

i=l 

m— 1 



tr J2]=o( M i-o.j l i) In (a jA -M i )~(Mi-aj^^ 1 ) In (a^-i-Af;) 



n 

i=0 
m+1 

i=l i=l 

This model shares features with the Eguchi-Yang matrix model [55], see also |13j . 



i=0 

m+1 m 
i=l i=l 



6 Conclusion 

We have rewritten the topological vertex formula for the partition function of the topo- 
logical A-model as a matrix integral. 

Having expressed the topological string in terms of a matrix model, we can bring the 
immense matrix model toolkit which has been developed since the introduction of random 
matrices by Wigner in 1951 to bear on questions concerning the topological string and 
Gromov-Witten invariants. We already started down this path in section [5] above. Going 
further, we can apply the method of bi-orthogonal polynomials [33] to our matrix model to 
unearth the integrable system structure (Miwa-Jimbo [211 EZ]) underlying the topological 
string, at least in the case of toric targets, together with its Lax pair, its Hirota equations 
(which arise as orthogonality relations), etc. In a related vein, free fermions [5H1 EH] 
arise in the theory of matrix models when invoking determinantal formulae to express 
the matrix model measure [60J. It will be very interesting to explore how this is related 
to the occurrence of free fermions in topological string theory, as studied in [HD E21 E31 
IM] . More generally, one should study what can be learned about the non-perturbative 
topological string from its perturbative reformulation as a matrix model, as in the works 
[65] 1661 167] 168] . A recurrent such question, which could be addressed in the matrix model 
framework (in fact, it was already latently present in the calculations in this work), is 
that of the quantization of Kahler parameters. 

On a different note, notice that the matrix model derived in this article, with a potential 
which is a sum of logs of g-deformed V functions, looks very similar to the matrix model 
counting plane partitions introduced in [69]. This is a hint that it could be possible to 
recover the topological vertex formula, corresponding to the topological string with target 
C 3 and appropriate boundary conditions, directly from the matrix model approach. Either 
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along these lines or the lines pursued in this paper, it would be interesting to derive a 
matrix model related to the Nekrasov deformation [TTJ |7U] of the topological string. 

A completely open question is whether the close relation between topological strings and 
matrix models persists beyond toric target spaces, and more ambitiously yet, whether 
there exists a general notion of geometry underlying matrix models. 
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A q-product 

The g-function, which plays a central role in the definition of our matrix model, is defined 
as an infinite product, 

g(x) = H(l--f). 

n=l 

It is the quantum Pochhammer symbol g(x) = [q/x;q]oo, and it is related to the q- 
deformed gamma function via T q (x) = (1 — q) l ~ x g(l)/g(q 1 ^ x ). 

The RHS is convergent for \q\ < 1 and arbitrary complex x ^ 0. g{x) satisfies the 
functional equation 

g(qx) = (1 - -)g(x) . 
x 

For n G N, we have 

g(q n ) = o 

and 



n ~ 1 /1\2 
-i n(n+l) t — r n(n+l) , n(n+l) Qi 1 

g'{q n ) = (-1)^(1) g-V* Y[[m]=g(l)q-^ A [n-iy^i-ir-'q-^ J[ 



g(q l ~ n ) ' 
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Via the triple product representation of the theta function, 

oo 

0^. T j _ J^J(1 — e 2nimT )(l + e ( 2m - 1 ) 7Tir + 27Tiz ^l _|_ e C2m-l)niT-2niz^ ^ 
m=l 

we obtain the identity 



1 1,9 ln ?\ / x ,Qs 

2 + 4^ ln ^ ; 2^r 9{x)9{ x' 



We have 



g(x)g(q/x) e^) 2 -^ln^ — In q / In x 1 Z7T — 2i7r 

g 2 In q g In q — . . jj ' 



g(l) 2 Vx~ Vlng 2 lng' lng 

where 6 is the Riemann theta-function for the torus of modulus — 2iir/liiq. 
At small lng, the following expansion is valid, 

^9{x) = Y~ E (In?)" Li 2 _ n (l/x) , 

mo n! 

^ n=0 

where we have used the definition of the Bernoulli numbers B n as the coefficients in the 
expansion of tj (e* — 1), 

t °° t n 

e { — 1 n! 

n=0 

Li n is the polylogarithm function, defined as 

LUx) = E ^ ■ 

k=l 

This is a generalization of the logarithm function, recovered at n — 1, 

Lii(x) = -ln(l -x) . 

It satisfies the functional relation 

Li / n (x) = iLi I1 _ 1 (x). (A.l) 

.X. 

Note in particular that this implies that Li n is an algebraic function of x for n < 0. 
E.g., 

Lio(x) = ^ . 

We also define the function 



OO K 

X K 



g(x) 
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Using the functional equation (A.l) of the polylogarithm, we find its small ln(g) expan- 
sion 



1 oo 

— Y 



■l) n B. 



n=0 



111 



n Qnq) n Li x _ n (l/x 
lng 



B 



2n 



lno I 111 *- 1 2(.r - 1) ^ (2//)! 



(lng) 2n Li 1 _ 2n (x) 



For the second equality, we have used B = l,Bi = — |, and -E?2n+i = for n > 1. 
We have near x — > oo 



and near x — > 0: 
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